Origin of room-temperature ferromagnetism in Mn doped semiconducting CdGeP2 
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CdGeP2 chalcopyrites doped with Mn have been recently found to exhibit room temperature ferro- 
magnetism. Isovalent substitution of the Cd site is expected, however, to create antiferromagnetism, 
in analogy with the well-known CdTe:Mn (d*) case. However, chalcopyrite semiconductors exhibit 
low-energy intrinsic defects. We show theoretically how ferromagnetism results from the interaction 
of Mn with hole-producing intrinsic defects. 



Spintronics - the use of spin rather than charge in 
electronics - cal s for room-temperature ferromagnetism 
in a semiconductor. While there have been several at- 
tempts to dope semiconductors with transition metal 
atoms 0, most of them did not show magnetic order- 
ing. The dopant of choice here ||] is Mn° (d^s^), being 
formally Mn^+ (d^) when substituting a divalent cation, 
and Mn'^+ (d*) when replacing a trivalent cation. Mn 
on the divalent cation site of II- VI semiconductors is at- 
tractive because of the significant solubility of the iso- 
valent Mn^+, but this produces antiferromagnetism [Q. 
The reason is that in a ferromagnetic arrangement all 3d 
spin-up orbitals are occupied so there is no empty or- 
bital on a neighboring Mn atom for the electron from 
one Mn atom to hop onto. Therefore an antiferromag- 
netic arrangement of spins at neighboring Mn sites is fa- 
vored. The III-V compounds doped with Mn would have 
been suitable candidates for ferromagnetism as these sys- 
tems are known |^ to have holes (d^ d^ + hole) when 
Mn'^+ substitutes the trivalent cation site. However, the 
ferromagnetic transition temperature (Tc) in these com- 
pounds is limited by the fact that Mn'^^ has low solubility 
and above a critical concentration 8% in GaAs) 
tends to phase separate. This was not a limitation while 
doping the II- VI semiconductors. Recently, Medvedkin 
et al. came up with the novel idea of obtaining both 
high solubility and high Tc in a semiconductor: replacing 
Cd sites in the chalcopyrite CdGeP2 with Mn atoms. Al- 
though CdGeP2 is isovalent to the III-V compounds 
(the average valence of Cd and Ge is a group III element), 
a large amount 50% at the surface) of Mn could be 
doped into this system, resulting in ferromagnetism in a 
semiconductor at the unprecedented high temperature 
(Tc=320 K) 1^. However, the existence of ferromag- 
netism in CdGeP2:Mn is surprising since, as explained 
above, replacement of a divalent site by Mn^+ (d^) is 
known to give rise to an antiferromagnetic interac- 
tion between the Mn atoms. Indeed, recent electronic 
structure calculations |^ found that Mn doped CdGeP2 
was antiferromagnetic, just like Mn doped CdTe We 
show here that Mn doping in CdGeP2 can produce stable 
ferromagnetism. This hitherto unexplored unusual form 
of ferromagnetism results from the interaction of a mag- 
netic ion such as Mn with holes produced by intrinsic 
defects. The central point is that chalcopyrite semicon- 



ductors are known pO| to be stabilized by certain intrin- 
sic defects such as cation (Cd,Ge) vacancies, vacancy- 
antisite pairs and the presence of hole-producing defects 
which could result in ferromagnetism being favored even 
when Mn dopes the Cd site. The theoretical challenge 
is to identify the hole-producing defects that form stable 
complexes with substitution. In this work we have cal- 
culated the formation energies of various kinds of defects 
and predicted the conditions favoring the substitution of 
Cd sites and Ge sites by Mn. From this we infer that fer- 
romagnetism is induced by holes generated in structures 
resulting from the simultaneous substitution of both Cd 
and Ge sites or the substitution of only the Ge sites with 
Mn. This offers a novel design principle for obtaining 
both high Mn solubility and high- Tc ferromagnetism due 
to defect-induced hole production in semiconducting sys- 
tems. 

We first calculated the formation energies of intrinsic 
point defects such as a vacancy at a Cd site (Vcd), a va- 
cancy at a Ge site (Vce) and a Ge antisite on Cd (Gecd), 
finding the energetically stablest defects. Then we cal- 
culate formation energies for Mn substituting a Cd site 
(Mncd) and a Ge site {Mnce)- Finally we deduce which 
of the many possible complexes among such single point 
defects and Mn would have simultaneously low formation 
energy (hence, high concentration) and create hole states 
(therefore coupling the Mn atoms ferromagnetically). 

We performed plane-wave pseudopotential |pl]-p^ cal- 
culations with isolated defects/Mn atoms introduced into 
64 atom supercells of CdGeP2. The lattice parameters of 
these supercells were fixed at the theoretically calculated 
values for CdGeP2: a=5.81 A and c=10.96 A ; while the 
atomic positions were relaxed. The formation energy for 
a defect comprising of atoms a in the charge state q was 
computed using the expression |10| 



AHf{a,q) = E{a) - E{0) 
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where E{a) and E{Q) are the total energies of the su- 
percell with and without defect a. Ua is the number of 
each defect atom; Ua = -1 if an atom is added, while 
= 1 if an atom is removed. /i„ is the absolute value 
of the chemical potential of atom a. Since the forma- 
tion energies are conventionally defined with respect to 
the elemental solid (s), we express as the sum of a 
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component due to the element in its most commonly oc- 
curring structure /z^, and an excess chemical potential 
^ct *-e Ma = Ma + Ma- Here Ma for Ge, Mn and Cd 
are the total energies evaluated for the fully-optimised 
elemental solids in the observed crystal structures p4[ . 
If A H f{CdGeP2) is the formation energy of CdGcP2, 
then /icd and are determined by 

tiCd + f^Ge + 2np<AHf{CdGeP2) (2) 

Furthermore, fj,cd ^ 0; MGe < 0, because otherwise the 
elemental solids will precipitate. The presence of other 
intervening binary phases, however, further restricts the 
values of Mcd and ^ce- One must solve Eq. (2) along- 
with the constraints placed by the formation energies 
A Hf{Cd^P2) and A Hf{GeP) of Cd3P2 and GeP 

i^iCd + '2^lP<AHf{Cd3P2) (3) 
^iGe + ^Jip<AHf{GeP) (4) 

to find the allowed range for ^cd and MGe in CdGeP2. 
The electrons ionized upon forming a positively-charged 
defect join the Fermi sea so the formation energy in- 
creases by (7 £/ , where £/ is the fermi energy which varies 
from eV at the valence band maximum (VBM) of the 
host material to the band gap of the host [|l5| . Equations 
(2)-(4) were solved using the experimental values [ p^JT?! 
of the formation energies for the binary phases Cd3P2 
(-1.2 eV) and GeP (-0.3 eV), while a value of -1.5 eV, in 
the same range as other chalcopyrites |l^ was used for 
CdGePa. 

The allowed range of chemical potentials Mcd and mge 
for CdGeP2 and the binaries Cd3P2 and GeP are given 
in Fig. 1. There are three distinct chemical potential 
domains where CdGeP2 can exist. For brevity we rep- 
resent each region by one characteristic point: point A 
- (Cd rich; Ge poor); point B - (Cd rich; Ge-rich); and 
point C - (Cd poor; Ge rich). Figure 2 shows the forma- 
tion energies of the intrinsic point defects Gecd, Vcd and 
Vce as well as substitutional defects Mnge and Mn^d 
at the chemical potentials A, B and C of Fig. 1 as a 
function of the Fermi energy. The vertical dashed line 
denotes the GGA gap which is underestimated with re- 
spect to the experimental 1.72 eV gap. Transition points 
between charge states are indicated by solid circles. The 
defects can form acceptor states generating holes in the 
system [e.g. (-/O) or (2-/-)] or donor states that generate 
electrons [e.g. {0/+)]- We see that: (i)When CdGeP2 is 
grown at point C (Cd-poor and Ge-rich), and to a lesser 
extent at point B (Cd-rich and Ge-rich), the favored sub- 
stitution of Mn is on the Cd site. This substitution leads 
to a neutral charge state (no holes), thus to AFM. (ii) 
However, for n-type conditions (E^? near CBM) at point 
B, and at all Ep values at point A (Cd-rich; Ge-poor), the 
stablest site for Mn substitution is the Ge site. This sub- 
stitution forms both single and double hole-producing ac- 
ceptors, which can promote FM. However, under n-type 
conditions at point B, the charge state q=-2 is favored, 



which has no holes and therefore promotes AFM. Analy- 
sis of the wavefunctions for different eigenvalues revealed 
that the Mn 3d levels are located 2-3 eV below the VBM 
and so one spin channel is fully filled (5 electrons). While 
formally Mn^+ would have a d"^ configuration, the exis- 
tence of a fully occupied d^ implies that two electrons 
are captured from the states near e/. The levels which 
are emptied correspond to the two acceptor state transi- 
tions taking place at 0.5 eV and 0.57 eV above FiypM- 

(iii) Cadmium vacancies with the charge state q=-2 are 
easily formed (in n-type conditions) at points A and C. 

(iv) Ge vacancy is stable for n-type conditions of point 
A, though the charge state which has the lowest energy 
in these conditions has no holes, can promote ferromag- 
netism. (v) Ge-on-Cd antisite has high formation energy, 
and would therefore not have appreciable concentration. 

Having identified the hole-producing centers that can 
yield ferromagnetism, we next examine the predicted 
solubilities of isolated Mn. Our calculated formation 
energies for CdGeP2:Mn and similar calculations for 
GaAs:Mn show consistently lower values (for the appro- 
priate chemical potentials) in the former case, predicting 
higher Mn solubility: the lowest formation energy of sub- 
stituting a Ga atom with Mn in GaAs is 0.5 eV (under 
Mn-rich, Ga-poor conditions). In contrast, even in the 
worst-case scenarios, we find a value of A Hf 0.6 eV 
for Miicd in CdGeP2 under Cd-rich, Ge-poor conditions. 
For MuGe we find A Hf ^ eV under Cd-poor, Ge-rich 
conditions. 

Having established low substitution energies for Mn 
in CdGeP2 we next see if Mn tends to cluster or not. 
We calculated the pairing energy, S, of a pair relative 
to infinitely separated Mn atoms in the host, given, by 
S = [E{2 Mn) - E{0 Mn)] - 2[E{1 Mn) - E{0 Mn)] 
where E(n Mn) is the total energy of the supercell with 
n Mn atoms. A negative pairing energy would imply 
that the two Mn atoms would prefer to pair rather than 
stay far apart. The pairing energy for nearest neigh- 
bor Mncd and Mnce pairs were found to be -1-577 meV 
and -1-513 meV respectively for the unit cell considered, 
implying no clustering. On the other hand, Mn substi- 
tution of GaAs had a negative pairing energy, equal to 
-250 meV for Mn occupying nearest neighbor Ga sites. 
The corresponding energy for clusters of three or four 
Mn occupying the vertices of a tetrahedron surrounding 
a single As in GaAs was even more negative, indicating 
that the Mn atoms showed a tendency of forming clusters 
when introduced into GaAs. The Mn atoms doped into 
CdGeP2 did not exhibit this tendency. 

Having established that the isolated defects Mncd, 
Mnce and Ycd have the lowest formation energy un- 
der certain experimental conditions and for some charge 
states the latter two are hole-producing (and therefore 
potentially ferromagnetism promoting), and that Mnge 
pairs repel each other rather than cluster we now study 
various combinations of these point defects in search of 
stable, hole producing complexes. Two Mn atoms were 
introduced alongwith various intrinsic defects at various 
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positions of a 16 as well as 32 atom super cells of CdGeP2 
corresponding to 50% and 25% Mn respectively. The 
energies of the ferromagnetic and antiferromagnetic ar- 
rangements was computed for the fully-optimised super 
cell {i.e both cell external parameters as well as cell- 
internal atomic positions were relaxed). The formation 
energies of the defect complexes in the 32 atom unit cell 
are plotted in Fig. 3 as a function of ficd in panel (a), 
and as a function of /ige in panel (b). The nature of 
the magnetic ground state (A=AFM and F=FM) has 
been indicated in parentheses for the defect complex with 
the lowest energy. The vertical arrows denote the values 
of chemical potential for which there is a crossover in 
the type of defect with the lowest energy. We see that: 
(i) Under Ge-rich conditions, where isolated Mn prefers 
the Cd site (Fig. 2c), a pair of Mn atoms also prefers 
the Cd site (Fig. 3a), leading to an AFM ground state. 
For 25% Mn the antiferromagnetic state is lower in en- 
ergy than the ferromagnetic state by 105, 19, 13 and 15 
meV/Mn for Mncd-Mncd separations of 1st, 2nd, 3rd 
and 4th neighbors respectively. For the 50% doping, the 
energy difference changes from 124 to 23 meV/ Mn. The 
energy difference of 19 meV and 23 meV for 25% and 
50% dopings for the second neighbor positions is in rea- 
sonable agreement with the difference of 21 meV and 
35 meV obtained by earlier calculations (ii) Adding 
a Cd vacancy to the Mncd-Mncd pair lowers the energy 
for the Cd-poor range of Fig. 3a, but leaves the system 
as a whole antiferromagnetic. (iii) Under Cd-rich condi- 
tions, where isolated Mn prefers the Ge site (Fig. 2a), a 
pair of Mn atoms also prefers the Ge site (Fig. 3b), ex- 
cept at the very Ge-rich end. For 25% Mn the FM state 
is lower in energy than the AFM state by 111, 67 and 82 
meV/Mn for Mnce-Mnce separations of 1st, 2nd and 3rd 
neighbors respectively. For the 50% Mn case, the ferro- 
magnetic state is still stable by 131 and 158 meV/Mn 
for 1st and 2nd neighbors. The fact that the energy 
difference between the ferromagnetic and the antiferro- 
magnetic states is quite large even when Mn atoms are 
at 3rd neighbor separation implies that the magnetic in- 
teractions are long-ranged in this case, unlike when Mn 
dopes the Cd site, (iv) The combination of Mncd-Mnce 
pairs is also ferromagnetic, being comparable in energy 
to when Mn replaces the Ge sites for certain chemical 
potentials of Fig. 3b. 

In pure CdGeP2 the valence electron configuration is 
a^tp, where p denotes the dominant character of the level. 
This reflects the fact that Ge contributes 1 electron to 
each of the 4 bonds with P, Cd contributes ^ electron and 
P contributes | electrons; the GeP4 tetrahedron trans- 
fers one electron to the CdP4 tetrahedron, through the 
bridging P. The electronic structure of a transition metal 
substituting a cation site in a semiconductor can be de- 
scribed by the hybridization between the orbitals of the 
transition metal ion and the anion dangling bonds formed 
from the cation vacancy in the host semiconductor [^|. 
Figure 4 shows this for Mnge in CdGeP2. The levels 



with <| character on Mn (left) hybridize with the lev- 
els with the same symmetry on the Ge- vacancy dangling 
bond (right), thus forming a bonding "crystal field res- 
onance" (CFR) tcFR-\ located deep in the valence band 
with dominant Mn character, and antibonding "dangling 
bond hybrid" (DBH) toBH] with dominantly p charac- 
ter. A similar hybridization of the down-spin levels re- 
sults in the down-spin counterparts of the CFR and DBH 
levels. The CFR level with e symmetry exists too. These 
SCFR, tcFR levels can hold 5 electrons. When Mn is 
substituted on the Ge site, 4 of its 7 valence electrons 
replace those of Ge, leaving 3 electrons to occupy the 
tcFR'\ ccFflT- However, since tjjBH is higher in energy 
than ecFR], two electrons move into ecFRh leading to 
i^CFR]^CFR^) consequently, two holes (empty circles in 
Fig. 4) are generated in the DBH levels. The intraatomic 
exchange splitting on the Mn atoms is large, ~ 0.6-0.7 eV 
Jigf , and what is unusual here is that it induces a large 
negative (0.8-1.0 eV) exchange splitting on the DBH lev- 
els which have dominantly P p character. The direction 
of the exchange splitting on the DBH (spin-down below 
spin-up) is opposite to the direction of exchange split- 
ting on the CFR levels, and this can be understood by 
simple tight-binding arguments. If the up-spin states on 
both the DBH and CFR are occupied, then the energy 
gain as a result of hybridization is small and comes from 
the single electron occupying the down spin DBH levels. 
However, if the up spin states on the CFR and the down 
spin states on the DBH are occupied, the energy gain 
as a result of hybridization is huge and there are contri- 
butions from both up and down spin channels. Hence 
the configuration in which the up spin levels on the CFR 
and the down spin levels on the DBH are occupied is 
the ground state of the system. A similar hybridization- 
induced antiferromagnetic coupling was earlier proposed 
to explain the high Curie temperature observed in an- 
other ferromagnet in which the magnetic Fe atoms are 
separated by nonmagnetic Mo atoms ||2^. In the inset 
of Fig. 4 we show the contributions to the total DOS 
from the up-spin states as well as the Mn d states in a 
small energy window near the fermi energy. The states 
are found to be strongly spin-polarized, with just 25% 
Mn d character from the two Mn atoms present in the 
32 atom super cell considered here. As a result of the 
large exchange splitting induced on the DBH states, the 
downspin DBH bands are energetically degenerate with 
the dominantly P p bands. The strongly polarised as well 
as delocalised down spin DBH band we believe is respon- 
sible for the long-ranged magnetic interactions, and, con- 
sequently the room temperature ferromagnetism found in 
Mn doped CdGePs. 
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FIG. 3. The formation energies as a function of (a) Cd (b) 
Ge chemical potentials calculated for complex defects in a 32 
atom cell of CdGeP2. The magnetic ground state (A=AFM, 
F=FM) for the lowest energy state is indicated in parentheses. 
The vertical arrows denote the value of chemical potential for 
which there is a crossover in the lowest energy defect type. 

FIG. 4. The energy level diagram of Mnge in CdGeP2. 
The energy of the levels are indicated in brackets. The inset 
shows the contributions to the total DOS in the energy inter- 
val -0.5 eV to 0.5 eV from the up spin states as well as Mn 
3d calculated for 2 Mnce in a 32 atom unit cell of CdGeP2. 



FIG. 1. The range of Cd and Ge chemical potentials where 
CdGeP2, GeP and Cd3P2 are stable. 



FIG. 2. The formation energies of different charge states 
of isolated defects calculated as a function of fermi energy in 
a 64 atom cell of CdGeP2 for three chemical potentials: A, B 
and C. The vertical dashed line denotes the calculated GGA 
band gap. Solid circles denote transition energies between 
charged states. 
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